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ABSTRACT 

We constrain the mean kinetic efficiency of radio-loud active galactic nuclei by using an optically selected 
sample for which both the optical and the radio luminosity functions (LFs) have been determined; the former 
traces the bolometric luminosity L, while the latter traces the kinetic power L^, empirically correlated to the radio 
emission. Thus in terms of the ratio = LyjL, we can convert the optical LF of the sample into a radio one. 
This computed LF is shown to match the directly observed LF for the same sample if = 0.10!j]:|]^ holds, with 
a scatter a = 0.38!;;;jj^ dex; with these values we also match a number of independent correlations between Lk, 
and radio emission, that we derive through Monte Carlo simulations. We proceed to translate the value of gy into 
a constraint on the kinetic efficiency for the production of radio jets or winds, namely, ty = L^^/M^c^ ~ 0.01 in 
terms of the rate Ma of mass accretion onto the central black hole. Then, on assuming that on average the radio 
sources share the same kinetic efficiency, we compute a sohd lower limit of about 25% on the contribution of 
radio sources to the local black hole mass density. 

Subject headings: cosmology: theory - black hole: accretion - quasars: general - radio sources: evolution 



1. INTRODUCTION 

The origin of the radio-loudness of Active Galactic Nuclei 
(AGN) still constitutes an open issue. There is evidence sug- 
gesting that the formation of a relativistic jet or a fast wind 
(e.g., Blundell & Kuncic 2007) sustaining the radio emission 
is tightly related to the mass of the central black hole (BH; 
e.g., Laor 2000). On the other hand, the wide scatter observed 
between radio and optical luminosities (e.g., Cirasuolo et al. 
2003a) suggests that other parameters such as the mass accre- 
tion rate onto the BH and possibly its spin could play a sig- 
nificant role in determining when a galaxy becomes radio-loud 
(e.g., Blandford et al. 1990, Ho 2002, Best et al. 2005, Sikora 
et al. 2005 and references therein). 

At one enticing extreme, Blandford & Znajek (1977) pro- 
posed the jets to be powered by the extraction of energy al- 
ready accumulated in a rotating BH. On the other hand, the 
spin, however important to set the jet direction, may not provide 
the necessary power for energizing the very luminous sources 
such as some Steep Spectrum Radio sources (SSRs) and many 
Flat Spectrum Quasars (FSQs). 

Alternative models (e.g., Livio et al. 1999) have proposed 
the dominant fraction of the jet or wind kinetic power to be di- 
rectly linked to the rest-mass energy of the currently accreting 
matter, as is the case for the radiative power. In any event, the 
kinetic power must originate at some time from the accretion 
energy; thus a constraint on the link between these two quan- 
tities is of key importance to probe the still unsettled issues of 
the jet's composition and production mechanisms (e.g., Celotti 
2004). 

In this work we propose a simple but efficient and statistically 
significant way to constrain the fraction of kinetic to accretion 
power. In the case of the radiative bolometric luminosity L, the 
emitted power is written as L = eMaC^, where Ma is the current 
accretion rate of matter onto the central massive BH and e is the 
efficiency to extract L from Ma- In an analogous way we can set 



the kinetic output as Lk = e^kM^c^ = {£i^/e)L. We assume that the 
AGN kinetic output represents a fraction gkof the bolometric 
luminosity, i.e., 

Lk = gkL; (1) 

this implies that gk = €k/e if the two outputs draw from the same 
accretion flow. Thus one could constrain g^ and in turn by 
selecting a statistically significant sample of radio-loud AGNs 
for which both the radio and optical luminosities are known. 
The latter is in fact a good tracer of the bolometric emission, 
given that AGN spectral studies have provided over the years 
a reliable average bolometric correction Cb = L/LbVb ~ 10, 
where the B-band luminosity Lb is in units of ergs~'Hz"' and 
vb = 6.8 X 10''* Hz (Elvis et al. 1994). In turn, significant radio 
emission is always associated with kinetic power, and several 
papers (referred to in § |2]i provide and discuss empirical corre- 
lations between Lk and the radio luminosity Lr. 

Our aim is to statistically derive the fraction of the accreting 
rest mass energy that powers the jets or winds by computing the 
average ratio g^, and then derive a lower limit to the predicted 
local BH mass density associated with radio-loud AGNs. In §|2] 
we will present our method, the results of which will be given 
in § [3] In § |4] we will describe in detail how to compute the 
corresponding contribution to the BH mass density using the 
value of Ek ■ In § E] we will discuss our results and give our 
conclusions. 

In the following we will adopt the "concordance cosmol- 
ogy" (see Spergel et al. 2007) with round parameters Ho = 
70kms-iMpc-i, Vl„ = 0.3 and Ha = 0.7. 

2. THE METHOD 



In this §we describe how we statistically constrain the ratio 
between kinetic and bolometric luminosities gk- Our database 
will be constituted by the optical and the radio luminosity 
functions (LFs) derived by Cirasuolo, Magliocchetti & Celotti 
(2005; henceforth CMC) for the same sample of optically se- 
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lected radio-loud AGNs; this is collected by cross-correlating 
the 2-degree Field Quasar Redshift survey (2dF) with the Faint 
Images Radio Sky at Twenty cm (FIRST). The sample spans a 
considerable range in magnitude -24 < Mg < -28 and redshift 
0.8< z ^ 2.2. In the following we will start from their optical 
Lpjj with parameters adjusted to our cosmology. 

We then convert the adopted optical LF to a bolometric one 
by using the average value Cb = 10.4 with a log-normal scatter 
of 0. 1 dex around the mean. We have checked that our results 
are not sensitive to the precise value of Cb, as very similar con- 
clusions are found also on using Cb « 8 (e.g., Marconi et al. 
2004, Richards et al. 2006). The direct link between optical 
and radio luminosities has been already extensively studied by 
CMC (and references therein) who find that a correlation be- 
tween these two quantities exists, although with a wide scatter. 

In this work we are mainly focused on constraining the av- 
erage kinetic output of radio-loud AGNs. Several studies have 
been able to empirically define the relation between radio emis- 
sions and kinetic outputs, showing that on average the former 
are associated to high levels of the latter (e.g., Rawling and 
Saunders 1991). After expressing the bolometric LF in terms 
of kinetic luminosities via Eq. ([T), we convert to radio powers 
on using the relation of Willott et al. (1999) 



U = 3x 10'' f^^L%] ergs- 



(2) 



is the observed radio luminosity in units of 
1 at 151 MHz. 



Our computed radio LF, shown as a solid line in Fig.[T] fits 
the CMC data on the empirical radio LF (solid points) when 
„ _ n 1 n+0 05 

a = 0.38to9- (3) 
The grey area in Fig. [T] shows the propagated la uncertainty 
from the optical LF. The match is good at all luminosities and 
redshifts in the sample using the same value for (shown are 
the results for z = 1 .3 and z = 2). 

We note that on adopting a Gaussian scatter we are able to 
reproduce the mean and extent of the observed scatter around 
it, within a given bin of optical luminosity. Fig. |2] shows our 
results (solid lines) for -24 <Mb< -25 and -25 <Mb< -26 
compared to the data collected by CMC (dashed lines). Ne- 
glecting scatter in the relations (dotted lines) would instead lead 
to a severe underestimate of the significant spread observed in 
the data. 

We have also performed other checks and tests on our re- 
sults. To look for possible biases in the optical selection of 
radio sources, we have compared our determination of the op- 
tical LF with the recent results by Jiang et al. (2007). These 
authors have cross-correlated a large sample of optically se- 
lected quasars from SDSS with FIRST, determining the radio- 
loud fraction of quasars as a function of redshift and luminosity. 
By multiplying the LF by Richards et al. (2005), representative 
of the whole population of optical AGNs, by the Jiang et al. 
(2007) fraction of radio sources as a function of optical lumi- 
nosity, we can evaluate the radio-optical LF expected from the 
SDSS data; in Fig. [3] we show this with a dashed line in two 
reference redshift bins at z = 1.3 and z = 2. Also shown with 
a striped area is the I- a uncertainty region derived from the 
uncertainty in the Jiang et al. (2007) radio-loud fraction. On 
the same plots we compare the completeness-corrected evalu- 
ations by CMC for their radio-optical sample; we show their 
data points as solid filled circles with error bars, together with 
their best-fit estimate represented with a solid line. The good 
match between the two evaluations supports the absence of any 
specific bias in our adopted LFs, at least forM^ < -24.5 where 
most of the sample sources reside (cfr Fig.|2]i. 

Through Monte Carlo sampling we then extract a sample of 
sources from the optical LF with a given Lb 0, from which we 
compute L, Lk and Lr by using Eqs. ([T]i and (|2]i with their scat- 
ters having fixed ^k to the value given in Eq. (O. In this way 
we build simulated distributions and correlations among these 
observables, which we compare with a number of different data 
sets. 

Fig. |4](left panel) shows our simulated LR-Lk correlation, 
which is equivalent in slope, and slightly lower in normaliza- 
tion but still compatible, relative to the one derived by Heinz, 
Merloni & Schwab (2007), and shown by the long-dashed line 
with 1(7 uncertainties shown as dotted lines. Note that their re- 
lation is calibrated on a sample of FSQs, which could be more 
beamed and/or energetic than our SSR dominated sample (see 
also § |5]l. Also shown with solid points and error bars are the 
sets of data by Birzan et al. (2004) and Allen et al. (2006). 

As Eq. might appear to be model-dependent, in the right 
panel of Fig.|4]we also show how our results on Lk relate to the 
narrow emission lines luminosities LnlkH and compare with 
the data points collected by Rawlings & Saunders (1991). The 

^ Here we point out that the correct values of the CMC double power-law LF should read: a = 0.228, /3 = 2.37 for the slopes, $0 = 6.1 X 10"^ Mpc"^ for the 
normalization and M* (z) = M*(0)-2.5{kiz+k2Z^), with ki = 2.24, k2 = -0.696, for the break luminosity. 

* All simulation results presented in the following are provided by averages over 1000 realizations, each one with 10000 points. 

^ We use the conversions Lnlr = 3(3Loii + l.SLoill). = ^OIII and Lqii = Z./5 x lO"' with scatter equal to a (e.g Willott et al. 1999). 



where L151 
lO^'^WHz-i sr 

In Eq. (|2]l Lk was empirically measured on tracing the inte- 
grated pdV work done by radio AGNs in excavating the cav- 
ities observed in the hot gaseous medium around them. The 
relation calibrated by Willott et al. (1999) relies therefore on 
specific assumptions on how to link the age of the source to its 
kinetic power. Following Hardcastle et al. (2007), the factor 
/ is introduced in Eq. (|2| to account for systematic underesti- 
mates of the true jet power that this technique may induce; for 
example, the quantity / also includes the coupling efficiency 
between the AGN kinetic output and the surrounding medium. 
In the following we will use (and discuss in § |5]) the average 
value / = 15, as measured for a sample of Faranoff-Rayleigh 
I and II radio galaxies (Hardcastle et al. 2007 and references 
therein). 

Note that both and L151 in Eq. ^ are calibrated on sam- 
ples mainly composed of radio galaxies; however, if the uni- 
fication model is to hold on average for all radio sources (see 
Urry & Padovani 1995), AGNs and radio galaxies are similar 
physical systems only observed at different angles; so we will 
take Eq. (|2|i to hold also for our sample of radio-loud AGNs. 
Throughout we will adopt an average AGN radio spectral slope 
an = 0.7, typical of the bright Steep Spectrum radio sources 
sampled by CMC (see also references therein). A Gaussian 
scatter with variance cr is adopted around the mean of both 
Eqs. ([B and (|2]). 

Thus we obtain a radio LF which depends on only two free 
parameters, gk and a, once / is fixed. The values of gk and a are 
then constrained through a analysis by matching our result 
to the radio LF independently measured by CMC. 

3. RESULTS OF OPTICAL-RADIO MAPPING 
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good agreement between our results and theirs, which were de- 
rived under different assumptions from Willott et al. (1999), is 
encouraging. 

4. THE KINETIC CONTRIBUTION TO THE BH MASS 
FUNCTION 

Constraining the total kinetic power from massive BHs is of 
key importance for several models of galaxy evolution (e.g., 
Croton et al. 2006, Cavaliere & Lapi 2007) which require a sig- 
nificant amount of kinetic feedback from the massive BHs to 
prevent catastrophic cooling in the cores of groups and clusters 
of galaxies and limit the formation of too massive galaxies in 
the local universe. 

Thus in this § we use our above results to evaluate the total 
kinetic energy associated to SSRs. As discussed later in § |5] 
we shall not include additional contributions to the kinetic in- 
tegrated AGN emission from other other kinds of AGNs, such 
as FSQs or radiatively inefficient sources, whose modeling is 
uncertain. Rather, our aim is to provide a firm lower limit to 
the actual contribution of radio-kinetic AGNs to the local mass 
density in black holes. 

To compute this, we reverse the classic Soltan (1982) ap- 
proach, which consists of estimating the average radiative ef- 
ficiency e by dividing the integrated energy density observed 
from AGNs by the local mass density in relic BHs. After com- 
puting the integrated kinetic energy from SSRs, we divide it 
by our best-fit value for the kinetic efficiency e^, and time- 
integrate to derive the local BH mass function expected from 
kinetically active SSRs. By comparing the latter with the local 
mass function of the whole BH population, we derive a lower 
limit to the fraction of the relic BHs which have been kinetically 
active in their past. 

We start from considering that during an accretion episode 
onto the central BH a fraction of the accreting mass energy is 
released as radiative and/or kinetic power In § [3] we evaluated 
the ratio gkto be of order 10"', which in turn implies a kinetic 
efficiency of 10"^, if the sources in the CMC sample pos- 
sess on average the same radiative efficiency e '--^ 0.1 typical 
of AGNs (e.g., Soltan 1982). Here we make the specific as- 
sumption that this value of is common to all Steep Spectrum 
kinetically-loud AGNs, irrespective of their radiative emission. 
The total energy extracted reads 

Ltot = L + Lk = (e + ek) MflC^ . (4) 

Therefore the total mass accreted onto the central BH is 

MBH = (l-e-ek)M„, (5) 

which converts to 



on using Eq. 0. 

By integrating the last equation over time and luminosity, and 
equating it to the local mass density pbh in relic BHs, we obtain 
the kinetic "Soltan-type" argument 

f dt f (l-e-ek)Lk 
J azj c^ek 
where $(Lk,z) is the kinetic luminosity function of AGNs in 
units ofMpc"^dex"' defined below. 

The value of derived in § |3] was calibrated on the CMC 
sample which is representative of radio-loud and optically 
bright AGNs, mostly composed of SSRs. Neglecting to low- 
est order e and ekin the numerator of Eq. (|7]i does not alter our 



conclusions. We now use as our wider database the 5 GHz ra- 
dio LP by De Zotti et al. (2005) specific for SSRs. The Steep 
Spectrum kinetic luminosity function <I>(Lk,z) is then derived 
by convolving this luminosity function with a Gaussian with 
mean given by Eq. ([T]i, where gk^nd the dispersion a are given 
in Eq. ©. 

Our result is shown by the solid line in Fig. |5] It is seen that 
at z « the BH mass density we evaluate to be associated to the 
Steep Spectrum radio-loud population is about 25% of the BH 
mass density found in local galaxies (light grey area). 

This amount is not far from the cumulative mass density ob- 
tained from the optical LF by Richards et al. (2005), represen- 
tative of the whole population of optical AGNs at z < 3, which 
is shown with a long-dashed line in Fig. |5] The latter mass 
density is obtained via the standard radiative Soltan argument 
which is formally equivalent to Eq. (Q with Ck in the denom- 
inator replaced by e, and Lk by L. Also shown with a dotted 
line in Fig.|5]is the BH mass density obtained by integrating the 
SSR LF by Dunlop & Peacock (1990); a slightly higher value 
would be obtained by integrating the SSR LF by Willott et al. 
(2001). These differences in BH mass densities simply reflect 
differences in the SSR LF adopted at given kinetic efficiency. 

In any event, our computations converge to show that the ac- 
creted mass density associated with Steep Spectrum radio-loud 
kinetic AGNs alone is comparable to the total BH mass den- 
sity accreted by optical AGNs. Note that the optical luminosity 
function has been extended down to Mb ~ -2 1 , which is the 
minimum luminosity probed by Richards et al. (2005). On the 
other hand, our mass density evaluations derived from the SSR 
LFs have been computed by extending the computations to the 
corresponding limiting luminosity, i.e. logLk/(ergs"') ~43.9 
and logP/(WHz"' sr"') - 23.1, on the basis of Eqs.|2]and [3] 

In the left panel of Fig. |6]we show the differential mass den- 
sity as a function of z for radio-loud and optical AGNs. It can 
be seen that the approximate match between the mass density 
accreted by kinetic and optical AGNs holds at each redshift, 
supporting a scenario in which the optical and kinetic outputs 
are produced by similar accretion events at each time. 

The match also holds in the final cumulative BH mass func- 
tions, as shown in the right panel of Fig. |6] The grey area in the 
Figure represents the local BH mass function with its system- 
atic uncertainties. Keeping in mind that below a BH mass M, ^ 
10** M0 our results are actually extrapolated, it can be seen that 
the fraction of the local BH population which has been kinet- 
ically active as Steep Spectrum sources in its past grows with 
M, from about 25% to reach near totality for masses > 10^ M©. 
This was computed on adopting an average Eddington ratio of 
m = L/LEdd = 0.5, where LEdd= 1-3 x lO^^^M./Mocrgs"'. This 
ratio is found to provide a good match between the local and ac- 
creted mass functions for the overall AGN population by, e.g., 
Marconi et al. (2004), Shankar et al. (2004) and Shankar, Wein- 
berg & Miralda-Escude (2007); our results however do not sig- 
nificantly depend on the precise choice for the adopted Edding- 
ton ratio, within the range of several tenths. 

The results presented in this § indicate a considerable con- 
tribution to the local mass density from radio-kinetically ac- 
tive SSRs, in fact comparable to the one from optical/X-ray 
sources (e.g., Shankar et al. 2004). The latter ones, not ex- 
plicitly accounted for here, are able to explain the mass around 
M, ^ 10** M0 needed to attain full match with the local mass 
function (see right panel of Fig. |6]l. In other words, a consid- 
erable fraction of relic black holes have possibly undergone a 
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kinetically-loud phase in their past. 

5. DISCUSSION AND CONCLUSIONS 

We have proposed a powerful method for estimating the aver- 
age efficiency of the kinetic power production process, derived 
essentially from observational constraints, i.e., from intercali- 
brating the radio and the optical LFs. Our result therefore is in- 
dependent of specific assumptions on the jet or wind dynamics 
and composition; it can be actually used to constrain theoretical 
models of jets. The value for the constant / ~ 15 discussed by 
Hardcastle et al. (2007) and inserted in Eq. (|2|i, supports the 
picture of "heavy" jets with a dominant protonic component; 
this enhances the kinetic output associated with even modest 
radio luminosities. We have adopted Eq. (|2|i for all radio-loud 
sources, irrespective of their redshift and/or environment. 

So far we did not consider effects of beaming in the Steep 
Spectrum radio LF. We expect the degree of beaming in these 
sources to be small on average. In fact, we have checked that 
on correcting the observed luminosities in the CMC radio LF 
with the beaming parameters appropriate for SSRs as done in 
the unification models of Urry & Padovani (1991, 1995), the re- 
sulting intrinsic LF is very similar to the one observed by CMC, 
only sa 0.06 dex fainter on average. This in turn implies values 
about 10% smaller for the best value, still within our I - a 
uncertainties. 

Our statistical results on the average kinetic efficiency of 
SSRs do not exclude that FSQs could have a higher kinetic ef- 
ficiency, or that there could be sources with exceptionally high 
kinetic outputs. For example, on discussing individual powerful 
blazars Maraschi & Tavecchio (2003) find higher values for the 
kinetic efficiency than derived here. This is partly due to their 
selection of the brightest jets for which the SEDs can be deter- 
mined up to 7-ray range. In addition, the high energy emissions 
probe the jet closer to the nucleus; conceivably the jet deceler- 
ates before reaching the radio hot spots and lobes, so that part 
of its kinetic energy is lost on the way. Thus the extended radio 
emissions on which Eq. (O has been calibrated could only set 
a lower limit to the full kinetic energy imparted to the jet at its 
origin. Heinz et al. (2007) have recently estimated the cumu- 
lative ESQ kinetic energy, finding average kinetic efficiencies 
a factor of a few higher than our best-fit value, but still on the 
order of a few percent of accreted rest-mass energy. 

To make contact with these independent studies, we have re- 
done our calculation on adopting for the correction factor / a 
higher value than the already significant one used here; for ex- 



ample, setting / = 45 in Eq. dU, we get a relation very similar to 
that derived by Heinz et al. (2007). Our fitting yields in this 
case gk=0.53![;:^^' and a = 0.307!|j {Jg^, which implies kinetic ef- 
ficiencies a factor of five higher than derived in §|3] While more 
stringent calibrations on relations such as that in Eq. (O are es- 
sential to pin down the average value of , all methods support 
kinetic efficiencies of order a few percent. Consider that having 
increased / and gk will proportionally increase the kinetic ef- 
ficiency and mean kinetic luminosity, although leaving the BH 
mass density pbh Lk/c/t unaltered. 

Finally, we note that the levels of kinetic efficiency derived in 
this work agree with the amount of energy kinetic feedback re- 
quired in theoretical studies of massive galaxies (e.g., Granato 
et al. 2004, Cavaliere & Lapi 2006, Shankar et al. 2006). 

Summarizing, we have used an optically selected sample of 
radio loud AGNs dominated by SSRs to convert the optical LF 
to a radio one. We have found that the kinetic output in jets or 
winds must amount to a fraction 10"' of the bolometric lu- 
minosity to match the empirical radio LF for the same sample; 
with this value we have been able to reproduce also a number 
of empirical correlations relating Lb, Lk and Lr. From typi- 
cal values for the AGN radiative efficiency e ^ 10"', we have 
derived an average kinetic efficiency e\^^ 10"^ which we have 
assumed to be common to all SSRs. By using Eq. we have 
then constrained the contribution of kinetically-loud SSRs to 
the local BH mass density to be > 25%, comparable to what is 
found from optical/X-ray selected AGNs. We consider the lat- 
ter estimate to provide a firm lower limit; additional contribu- 
tions to the kinetic integrated AGN emission could come from 
the more numerous and less luminous AGNs (e.g., Heinz et al. 
2007), from AGNs in radiatively inefficient but kinetically effi- 
cient accretion states (e.g., Churazov et al. 2005), and/or from 
FSQs with intrinsically higher kinetic power than here consid- 
ered (see Maraschi & Tavecchio 2003). An additional contribu- 
tion arises from radio galaxies which should have comparable 
kinetic outputs according to Unification models; in agreement 
with our findings in fact, Koerding et al. (2007) find gk ~ 0.2 
for a significant sample of FRI and FRII radio galaxies. 
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Fig . 1 . — Best-fit radio LF (solid line) and its propagated 1-a uncertainty region (grey area) obtained by converting the optical LF from CMC 
with the use of Eqs. ([TJ and The result is presented for two redshifts bins as labeled, and is compared with the data on the radio LF empirically 
derived by CMC (solid circles) for the same sample. 
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Fig. 2. — Comparisons between the normalized distributions of sources as a function of radio power corresponding to a given bin in optical 
luminosity, as labeled. The dashed line in each panel corresponds to the sample by CMC, while the solid line is the mean computed from our 
simulations with scatter included. The dotted line is the distribution expected when no scatter is considered. 
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Fig. 3. — The solid line (with its l-cr uncertainty region shown with a grey area) represents the fit to the optical LF, derived by CMC for their 
sample of optical-radio sources in the 2dF-FIRST catalogues; the solid points with error bars represent their original measurements. The dashed 
line and the striped area represent the best-fit and l-cr uncertainty region of the optical LF by Richards et al. (2005), representative of the whole 
optical AGN population, corrected for the luminosity and redshift-dependent fraction of radio-loud AGNs derived by Jiang et al. (2007). These 
two independent evaluations are compared at z = 1.2 and z = 2.0, as indicated in the two panels; note the good agreement. 
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Fig. 4. — Comparison between simulations and data. Left panel: the results from our best-fit model are compared with the set of data (filled 
circles) presented in Heinz et al. (2007) and with their model {Jong-dashed line and dotted lines). Right panel: our simulated best-fit relation 
(solid line with dispersion shown by the grey area) is compared with data derived by Rawling & Saunders (1991). All data have been converted to 
our adopted cosmology. 
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Fig. 5. — Cumulative black hole mass densities computed as a function of redshift. The long-dashed line represents the mass density as- 
sociated with optical AGNs, obtained by integrating the Richards et al. (2005)'s LF above Mb = -21 (i.e., logfi.4GHz/(W/Hz/sr) ~ 23.1 and 
logLk/(erg/s) ~ 43.9, on using eq. ([3}). The solid and dotted lines are the results of integrating above the same luminosity threshold the LFs by 
De Zotti et al. (2005) and Dunlop & Peacock (1990), respectively. The value of e = 0.1 and our best-fit value ^k=0.10 have been used for all 
curves. The boxes show the local black hole mass density pBH in massive spheroids (dark grey area) and 30% of this value (light grey area). 
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Fig. 6. — Left panel: differential accretion density as a function of redshift. Right panel: the corresponding cumulative number density at z ~ 
of black holes derived from the equation representing the building-up of p.(z) by gravitational interactions of their host galaxies (see Vittorini, 
Shankar & Cavaliere 2005, and Shankar et al. 2007). The line styles in both panels have the same meaning as in Fig. [5] The grey band shows the 
cumulative local mass function with its uncertainty region taken from Shankar et al. (2007). 



